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We study analytically the local density of states (LDOS) in a disordered normal-metal wire (N) 
at ballistic distance to a superconductor (S). Our calculation is based on a scattering-matrix ap- 
proach, which concerns for wave-function localisation in the normal metal, and extends beyond the 
conventional semiclassical theory based on Usadel and Eilenberger equations. We also analyse how 
a finite transparency of the NS interface modifies the spectral proximity effect and demonstrate that 
our results agree in the dirty diffusive limit with those obtained from the Usadel equation. 
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I. INTRODUCTION 

It is widely acknowledged that a piece of a normal 
metal that is in good contact with a superconductor 
acquires some superconducting properties. This phe- 
nomenon, named the proximity effect, has already been 
studied by Coopertl in the early sixties. Since then many 
theoretical and experimental investigations have been 
carried out.El Much owed to the recent progress in the 
fabrication technology of nanostructures there is a r&j 
vived interest to the proximity effect in the last decade.tl 
One remarkable evidence of this effect is the formation 
of a spectral gap in the normal metal, which strongly 
affects the low-temperature transport properties of the 
normal-metal superconductor (NS) junctions. The key 
mechanism responsible for the appearance of the gap is 
the Andreev reflection at the NS boundary, which con- 
verts the dissieative electrical current into dissipationless 
supercurrent.Q Similar mechanisms act in superconduc- 
tor ferromagnet junctionfcrwhich have become an object 
of intense study recently.QEl 

An effective experimental technique which allows for 
spatially resolved measurements of the electronic den- 
sity in the nanostructures is the scanning tunnelling mi- 
croscopy. It provides both a unique sub-meV energy sen- 
sitivity and an atomic spatial resolution. Several recent 
measurements of the local jX^lrjcjjxuiic density of states 
(LDOS) in the NS junctionsElffiBy turned out to be in 
very good flOTeernentjaqth the predictions of quasiclassi- 
cal theorjil3'E3oli3iljO of "non-equilibrium" supercon- 
ductivity,|_based on the Usadel equation for the diffusive 
transportlla^nd the Eilenberger equation for the ballistic 
transport 

The interplay of ballistic and diffusive transport be- 
comes important when one studies local properties at 
short distance to an NS interface in a disordered sys- 
tem. Quasiparticles are then transferred to the interface 
by ballistic transport, while they explore the rest of the 
system diffusively. This situation is not covered by con- 
ventional quasiclassical theory. Quasiclassics also cannot 
account for the non-perturbative effects of wave-function 
localisation, which only can be included by a fully phase- 
coherent approach. 

In this paper we present a theory that goes beyond the 



quasiclassical description and apply it to calculate the 
local density of states in an NS wire geometry near the 
interface, at zero temperature and vanishing magnetic 
field. 

In our model the normal metal is shaped in the form 
of the long disordered quantum wire, which supports N 
propagating modes at the Fermi level Ep. The elastic 
scattering mean free path i in the wire is assumed to be 
much larger than the Fermi wave length Xp, which cor- 
responds to the weak disorder. The superconductor is 
assumed to be clean and characterised by the bulk value 
A of the amplitude of the pair potential. The supercon- 
ductor order parameter is assumed to be constant A in 
the superconductor and zero in the normal metal. This 
approximation is refeped to in the literature as a "rigid 
boundary condition" .lj 

We calculate the mean LDOS, or, more precisely, its 
envelope, at a distance x on the normal-metal side of the 
NS interface as shown schematically in Fig. |l| The enve- 
lope is obtained by averaging the LDOS over distances of 
the order of the Fermi wave length Xp . The spatial aver- 
aging smears out the Friedel type oscillations and makes 
the LDOS independent on the position across the wire. 

We study in detail the case that the distance x is small 
compared to the scattering mean free path i, so that 
Xp ^ X <^ £, while the ratio between the superconduc- 
tor coherence length ^ — hvp/A and £ remains arbitrary. 
The resulting mean LDOS found by averaging over dis- 
order does not depend on x and is a smooth function 
of energy everywhere except at e = A (the energy e is 
measured from the Fermi surface). 

Our calculation is organised as follows. In Sec. H we 
derive a general relation between the one-point Green 
function in a quantum wire and the reflection matrices 
rp, rji. These matrices relate the plane- wave components 
of the quasiparticle wave function in the process of reflec- 
tion from the parts of the wire to the left and to the right 
part of X. 



We apply this result in Sec. Ill in order to calculate the 



mean LDOS in the neighbourhood of an ideally trans- 
mitting NS interface. The matrices r^, rp of the size 
2N X 2N describe the reflection of the electron- and hole- 
like quasiparticles. The left reflection matrix rz, is diag- 
onal in the electron-hole representation and depends on 
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FIG. 1: The geometry of an NS junction consisting of a 
long normal-metal disordered wire N, a clean superconductor 
S and a dielectric tunnel barrier I in between. The mean local 
density of states (LDOS) is calculated at the distance x from 
the NS interface, with Xp <^ x <^ £. The matrix vl relates 
the plane-wave components in the process of reflection from 
the normal-metal disordered wire. The matrix vr describes 
the reflection from the tunnel barrier - superconductor part 
of the junction. The mean LDOS is found by averaging over 
the disorder-induced fluctuations of the matrix vl- 



the disorder in the normal metal. The right reflection 
matrix rj^ is ofl^-diagonal (in absence of the tunnel bar- 
rier) and is fixed within the model considered. 

In the region e < A we obtain the disorder-averaged 
LDOS 



n{x,e) = T:pe{(j)A), 



arccose/A, 



(1) 



where the function ((/>) is the probability density of the 
eigenphase of the matrix correlator ro(£)ro(— e)^. The 
reflection matrix tq^e) relates the plane- wave amplitudes 
of the electron wave function in the process of reflection 
from the semi-infinite normal-metal wire. The probabil- 
ity density Pei4') has been studied in Ref. |2|. Apart from 
energy and phase it depends on the number of channels 
N and the mean scattering time Ts = £/vf- According to 
Ref. one can distinguish localised, diffusive and ballis- 
tic regimes in the form of the function Pei4>) depending 
on the value of e. We observe the effect of Anderson 
localisation in the linear increase of the LDOS for ener- 
gies smaller than the Thouless energy Sc = h/N'^Tg. We 
also find that the curves calculated for different number 
of channels in the wire are lying close to each other at 
any ratio (see Fig. ||). This suggests that the weak- 
localisation correction to the LDOS is small in the case 
of the ideally transmitting NS interface. 

In Sec. [V we generalise the model to include a tunnel 
barrier at the interface, parametrised by a tunnel proba- 
bility per mode F. We calculate analytically the LDOS 
near the interface in the extreme cases of a localised wire, 
N = 1, and a diffusive wire, N ^ 1. 

The effect of the tunnel barrier consists of a reduction 
of the pseudogap in the normal metal. This effect is most 
pronounced in the dirty regime ^ ^ or Ats/H ^ 1. 
The results of our calculation for the diffusive wire in 
the intermediate regime i — ^ are summarised in Fig. ^ 
for different values of F. We observe that the LDOS 
increases monotonously to its bulk constant value around 
the energy hT"^ /ts and reveals a high and narrow peak 



close to e = A. 

The monotonous reduction of the pseudogap is at- 
tributed to the quasiparticles which experience normal 
reflection at the tunnel barrier and therefore do not see 
the NS boundary. The formation of the peak is due to 
the quasiparticles reflected from the superconductor. 

When the distance x increases beyond the mean free 
path £ a competing effect takes place. That is the sup- 
pression of the pseudogap due to the back-scattering on 
the weak disorder potential in the normal-metal segment 
of length X in front of the interface. The estimated size 
of the pseudogap due to this effect is hD/x^, where D is 
the diffusion constant in the normal metal. In this case 
the LDOS considerably overshoots its bulk value around 
e = hD/x^, which is in contrast to the monotonous in- 
crease due to the tunnelling into the superconductor. We 
therefore anticipate that the effect of the tunnel barrier 
still can be seen in the shape of the LDOS provided 
hD/x^ » HT-^/ts, or equivalently x ^ i/T. Namely, 
at distances smaller than £/r the LDOS may acquire the 
step-like feature at the value e = HT'^/ts, which is fixed 
by the NS interface transparency rather than by the dis- 
tance to the interface. 

A qualitatively similar phenomenoB, has been indeed 
observed in experiments by Levi et. aln in the Cu barrier 
pin wires near a N(Cu)-S(NbTi) boundary. 

On the contrary, at large distances x » £/T the bar- 
rier is not effective in the sense that its presence cannot 
be distinguished in the energy dependence of the LDOS^ 
This is consistent with a general semiclassical criteriontj 
which states that the barrier is not effective for a given 
observable if the most of the relevant trajectories hit the 
NS interface more than F~^ times before the electron- 
hole coherence is lost. In the case of the LDOS this cri- 
terion is fulfilled ior X ^ £ /T . 

The tunnel barrier acts differently for the single- 
channel wire. In the dirty regime £ ^ ^ the size of the 
pseudogap HT/ts scales linearly with F due to the An- 
derson localisation. This results in a different shape of 
the LDOS compared to the diffusive case {N ^ 1). The 
difference becomes more and more pronounced with de- 
creasing ratio £/^ or tunnelling probability F. 

In Sec. we compare the LDOS for the diffusive case 
{N ^ 1) found from our tiieory to the LDOS calculated 
from the Usadel equation.llil 



II. GREEN FUNCTION IN A WIRE 
GEOMETRY 

In our model of the NS junction the normal metal is 
shaped in the form of a semi-infinite quasi-one dimen- 
sional disordered wire. The properties of such a sys- 
tem is welLunderstood in the framework of the scatter- 
ing theorjo provided the weak disorder limit Xp <^ £■ 
The detailed statistical description of the disorder scat- 
tering is based on. tha Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equationOnJ This is a scaling equation for the 
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probability distribution of the scattering matrix of a seg- 
ment of the wire. Below we derive a general relation 
between the one-point Green function and the reflection 
matrices r^, for two parts of the wire. The single- 
channel counterpart of this relation has been used re- 
cently to reconsider the^M-oblem of LDOS fluctuations in 
ID normal-metal wiresHj 

The disordered wire has the Hamiltonian H ~ 
Ho + V{r), where V{r) is a disordered potential. We 
parametrise f = {x,f}), where x is the coordinate along 
the wire and p is the vector in the transversal direction. 
We first discuss the case of 'spinless' electrons, assuming 
ffo = —2^^^) ?i = 1, and include hole- like quasiparti- 
cles in Sections III and IV. 

In the absence of V the quantisation in the transver- 
sal direction gives rise to a set of TV propagating modes 
characterised by the transverse momentum (f„. The 
total energy E = (l/2me)(|(fnp -I- fc^), where the x- 
momentum fc„ is conserved. The retarded Green func- 
tion G^{E) — {E + irj — H)^^ is written in the channel 
representation as 

G^„(x,x') = JJjpdp'{n\p){ff\m){r\G^\r'), (2) 

where the integration is carried out over a cross-sectional 
area A. Hence the LDOS 

n{f,e) = --Im V(m|p) (pl7i)G^,„(x, x), (3) 



where e is the energy measured from the Fermi sur- 
face. For a two-dimensional wire of the width d we 
have {p\n) = {2 / sin (nnp/d). In what follows we 
shall omit the index i?, assuming everywhere the retarded 
Green function. 

Let us formally cut the wire in the point x into two 
pieces and treat the left and the right part separately. 
We decompose the potential V — Vr + Vl, where Vb..l 
is the disorder potential in the right and the left part of 
the wire, respectively. We also introduce the left and the 
right Green function as G_r,l = {E + ir] — Hq — Vl,r)~^- 
According to Fisher and Lee, Ref. |2^, we have 



GL,B.;nm{x, x) 



1 



:{Snrn+rL,R;nrn{x)), (4) 



where w„ = fc„/me is the channel velocity and r^jf are 
the reflection matrices from the left and the right part of 
the wire, respectively. 

The Green functions obey Dyson equations which can 
be written in the matrix form as 



G{x,x) ^ Go{x,x) + dy Goix,y)V{y)G{y,x), (5a) 

^ — CXj 

G{x,x) ^Gr{x,x) + f dy GR{x,y)VLiy)G{y,x), (5b) 

J —OC 

poo 

G{x,x) ^ Gl{x,x) + dy GL{x,y)VR{y)G{y,x), (5c) 



where the elements of the matrix V are given by 

Vnm{x)^ f dp {n\p){p\m)V{r), 

J A 



(6) 



and the ballistic Green function (in absence of the poten- 
tial) reads 



Go,nm[X,X ) = e 



_ Onm ikn\x-x'\ 



(7) 



We also take advantage of the following relationsEEl 

GR,ni{x,y) = e"*(=""^)Gfl,„;(a;,a:), for y < (8a) 
GL,„i(x,y) = e'''''<^'^-y^GL.ni{x,x), fory>x, (8b) 

in the disorder-free regions in order to eliminate the in- 
tegral terms in Eq. (|^) . As a result we obtain the matrix 
equality 



1 



G(a;, x) 
Using Eqs. (E 



Go{x,x) Gr{x,x) 
and (|^) we finally get 



Gl{x,x) 



(9) 



G(x,x) = ^(1 + r;^)- 

V «V -L 



1 



rhrR 



-il + f^)— (10) 



where v is the diagonal matrix of channel velocities ?;„. 
Together with Eq. (H) this equation defines the LDOS 
via the reflection matrices. 

In the case of uncorrelated disorder the reflection ma- 
trices rij and rR are statistically independent, which 
makes Eq. (1^) useful for practical calculations. 

In genera the LDOS oscillates on the scale of Af (due 
to the prevailing contribution of one particular quantum 
state). These Friedel-type oscillations can play a crucial 
role especially in one dimension. In what follows we are 
concerned with the smoothed version of the LDOS that 
does not change on the scale of the Fermi wave length 
and, therefore, also not in the transversal direction. For 
this purpose we introduce the spatially averaged LDOS 



n{x, e) = 5V ^ n{r, e) dr, 
Jsv 



(11) 



where the integration is carried out over a volume SV 
around the point {x,p). The linear size of the volume 
SV is assumed to be much larger than the Fermi wave 
length and much smaller than the mean free path £. For 
— x'l <C ^ the reflection matrices defined at the cross 
section x' are related to those defined at x by 



rLix') 
rR{x) 



^ — ik{x — x') 



rL{x)e 



-ik{x — x') 



^ik{x — x') 



rR{x)e 



ik{x — x') 



(12a) 
(12b) 



with k — rngV. Expanding the right-hand side of Eq. 
( pO| ) in a geometric series in r^, tr we notice that only 
the terms with equal numbers of and rn matrices 
do not oscillate on the scale of the Fermi wave length 
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and have to be kept. Additionally the averaging in the 
transversal direction mixes up the different modes so that 
{m\p){p\n) cx Smn in Eq- (§)■ As the result we obtain 



n(x, e) = — ^ Rc Tr — 



(13) 



where no is the bulk value of the LDOS in the normal 
metal, which is set to unity in the rest of the paper. In 
what follows we apply Eq. (|l^) to calculate the LDOS 
in the normal-metal wire in the immediate vicinity of an 
NS interface. 



III. LDOS NEAR THE IDEAL NS INTERFACE 

The relation (|l^) applies straightforwardly to the 
model of the NS junction discussed in the Introduction. 
The only modification is the doubling of size of the re- 
flection matrices due to particle-hole conversion. We still 
denote the number of electron channels in the wire by N , 
so that the size of the particle-hole reflection matrix is 
now 2N . Equation (|l3|) can be written in the form 



2 

n{x,e) = 1 + — ReTr y^{rLrRY\ 



n=l 



(14) 



where is the electron-hole reflection matrix for the 
long normal-metal wire, while rn is that for the ideal 
NS interface. These reflection matrices are conveniently 
parametrised by 



rL 



roie) 
roi-ey 



1 

1 



(15) 



where (j)A — arccose/A is the Andreev phase and ro(e) 
[ro(e)*] is N X N reflection matrix of the electron- like 
[hole-like] quasiparticles for the normal-metal wire. The 
matrix product r^rR is block off-diagonal, hence only the 
even powers n contribute to the trace in Eq. (p^). From 
Eqs. (|l|,^ we obtain 



2 

n(x,e)-l + -ReTr^ 



(16) 



The right-hand side of Eq. ( p^ ) is completely deter- 
mined by the eigenvalues of the correlator ro(e)ro(— e)*, 
which is a unitary matrix. Its eigenvalues are conve- 
niently parametrised by exp(2j(^j), j = 1, 2, . . . , A^, where 
the phases are restricted to the interval (0, tt). The 
joint probability density Pi;{(j)i,(j)2, . . . ,4>n) is a symmet- 
ric function with respect to any permutation of its argu- 
ments because of the statistical equivalence of the chan- 
nels. This function has been studied in detail in Rcf. 

Our calculation is restricted to the mean LDOS 
n{x,s) = {n{x,s)), where the angular brackets corre- 
spond to the average over the disorder potential in the 
wire. In order to perform the average in Eq. (|T^), it is 
enough to know only the probability density Pe{<ti) of a 
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FIG. 2: The mean LDOS (|2^) in an N-channel normal-metal 
wire near an ideally transmitting NS interface. The curves 
are calculated from Eqs. ( |24| , p5| ) . The thick solid (dotted) line 
corresponds to the limiting case of the multi(single)-channel 
wire. The thin lines are for the finite number of channels 
A = 2, 3, and 4. The figures correspond to the clean regime 
Ats = 3 (top), the intermediate regime Ats = 1 (middle), 
and the dirty regime Ats = 0.3 (bottom). 



single eigenphase. It is instructive to compare Eq. ( p^ ) 
with the similar representation of the integrated density 
of states in the case of the normal-meta]-|Wire of finite 
length, which has been analysed recently .ca 

When the Andreev phase (pA is real, i.e., for e < A, 
the mean LDOS is found from Eq. (|l6|) as 



n(x,e) 



-KPe 



e < A, 



(17) 
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where the eigenphase density Pe(</') is assumed to be nor- 
mahsed to unity on the interval (0,7r). The probabil- 
ity density Pe{<t>) acquires its simplest, form in the case 
^ 1 of a large number of channelsEJ, 

pM = — iy- Im J(£T,)2+i£T,(l-e-2»^), (18) 
TT sm 6 * 



and in the single-channel cas 



TT 



exp {-£Tst) 



sin^ (f>-t sin 241+1 



-dt. 



(19) 



The scattering time of the DMPK scaling equation 
differs by a numerical factor (dependent on the dimen- 
sionality d of the Fermi surface) from the mean scattering 
time of the transport theory r^. Namely, Ts — CdT^^ where 
Cd = 2, 7r2/4, 8/3, for the dimensionality d = 1,2,3, cor- 
respondingly. 

Note, that the integrated density of states (DOS) v — 
L^^ n{x, e) dx in the iafinite disordered wire L ~* 00 
is given by the relationE3 v = n{d/de)[epg{0)], which 
is similar in spirit to Eq. (p^. For wires with on-site 
disorder (in standard universality classes) the value of 
7r/9e(0) equals to unity irrespective of energy e; however 
it can have a singularity at e = for wires with a specific 
disorder symmetry. 

So far we were only concerned with the mean LDOS 
for e < A. However, the result (17) can be easily ex- 
tended to the energies above the pair potential value 
with the help of the analytical continuation e — ioj. 
On the other hand the analytical continuation has an- 
other crucial advantage. It transforms the dynamical 
correlator ro(£)ro(— e)* into the essentially static ob- 
ject ro(it<;)ro(iw)*. In the absence of a magnetic field 
the time-reversal symmetry is preserved and the reflec- 
tion matrix tq is symmetric, hence ro(iw)* = ro(iCi;)^. 
The eigenvalues exp(2j0j) of the matrix ro(£)ro(— e)* 
are transformed to the real eigenvalues Rj of the matrix 
ro(ict;)ro(ici^)*, which are the probabilities of the reflec- 
tion from the long disordered wire in the presence of a 
spatially uniform fictitious absorption lo. 

The summation in Eq. ( p^ is performed in terms of 
the eigenvalues 



n(x, e) 



1 J^ l-flX(^) 

N ^ l + i?.a2(^) 

.7=1 ■' 



(20) 



-46-1-0+ 



where 0"*" is an infinitesimally small positive imaginary 
part of energy which ensures the retarded Green function 
required in Eq. (0). We have also introduced 



a{bj) 



^1 + {lo/A)^ -Lu/A. (21) 



The joint probability density of the eigenvalues Rj for 
the infinitely long wire is given by the stationary solution 
of the DMPK equation. In the parametrisation 

a, e (0,00), (22) 



this solution takes the simple form 

N 



k>j 



(23) 



which we recognise as the_prthogonal Laguerre ensemble 
of random matrix theorycJ (with normalisation constant 
Cm)- This ensemble corresponds to the class CI in the 
classification scheme of Ref. The probability density 
(one-point function) p(cr )_normalised to unity in the in- 
terval (0, cx)), is given byEa 

\n=0 

+ \L^^_,{a) j\ce^^-^y-L^^_,{C)y (24) 

where L^'if' [a) is the generalised Laguerre polynomial. 

We substitute the parametrisation ( p2| ) in Eq. (^0|) and 
average over disorder with the help of the density P({cr}). 
The result reads 



n(a;,e) — Ke da p{a) 



1 



2{N+1)ljt, 



(25) 



This equation extends Eq. ( |l7|) to energies larger than 
A. It can also be applied for arbitrary N. In the lar|ge-7V 
limit the distribution p(cr) can be approximated byEj 



lim NpiCN) ^—J--l. 



< C < 4. 



(26) 



Substituting this expression into Eq. ( |25| ) we reproduce 
the results of Eqs. (|l|,|l^) for e < A. In the limit e — > 
this leads to the square root behavior of the LDOS 
n(x, e ^ 0) — Re ^/—isTs. In the extremely dirty regime 
Ats — > we reproduce the result of the conventional BCS 
theory 



n{x, s) = Ree/ V 



A2 



(27) 



The Thouless energy Ec = 1/iV^Ts, however, remains 
unresolved within the multi-channel approximation (|2^) . 
In order to fix the scale fn one has to take advantage of 
another limiting relatione^ 

lim p(C/Af) = J2(2v/C)- Jo(2VC)J2(2VC) 

+ (2^y' Joi2^)M2^), (28) 

where Jn{z) are Bessel functions. In the limit e — > one 
can safely put a(w) = 1 in Eq. ( p5| ) and take the real 
part explicitly. 



CTj +2{N + 1)u;Ts 



00 

lix,s)=Jdap{a)^^-^^^ 



(29) 
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To leading order in e/ec the function p{a) in ( p9| ) can 
be approximated by its value at the origin p{0) — 1/2, 
which holds for a ^ N^^ [see Eq. (28)]. We therefore 
obtain 



n{x, e) 



7r(iV + 1) 



-£Ts 



^ I I 



), e<ec. (30) 



The factor in the last expression reflects the fact that 
only a single channel is responsible for the non- vanishing 
LDOS at energies lower than Sc- 

In Fig. I we plot the mean LDOS given by Eq. (H) 
against the ratio e/A for different numbers of channels in 
the moderately dirty regime Ar^ = 0.3, the intermediate 
regime At^ = 1, and the moderately clean regime Ar^ = 
3. We observe that the curves are lying close to each 
other in all cases. (This suggests that the LDOS near 
the ideally transmitting interface is quite insensitive to 
phase-coherent effects.) The situation changes in the case 
of a finite transparency F < 1 of the NS interface. 



IV. EFFECT OF A TUNNEL BARRIER 
A. Model 

We now introduce the simplest model of a dielectric 
tunnel barrier at the ideal NS interface. The mean LDOS 
is calculated in the normal-metal at a ballistic distance 
X from the interface (see Fig. |l|). 

We describe the segment / of the wire between the cho- 
sen cross section and the ideal NS interface (this segment 
includes the tunnel barrier) by its 5-matrix 



Si 



t 



(31) 



where each block itself consists of block-diagonal matrices 
in the particle-hole representation, 



1,2 



ri,2(-£)' 
ti,2(-e)* 



(32a) 
(32b) 



and the matrices ri.2(£), ^i,2(£) are N x N electron re- 
flection and transmission matrices corresponding to the 
segment /. 

The right matrix r/? in the fundamental formula ( |l^ ) 
depends on the S'-matrix of the segment / [see Eqs. 
(^,^2|)] and on the scattering matrix for Andreev reflec- 
~5|)]. A straightforward algebraic calcula- 



tion [see 
tion gives 



rR = 



trie) rrX-e)* 



(33a) 



tc{e) 
rc{£) 
M{e) 



e-''^'^('h2{-erM{e)h{e), (33b) 
n{e) + e-^''^^(^'h2{e)r2{-erM{e)h{e), (33c) 

l-e-2**^(^V2(£)r2(-e)*l"\ (33d) 



In general, if the segment / contains some weak disor- 
der (which is the case, for example, for x > £) the cor- 
relations between the matrices ri_2 and ti.2 for electron- 
and hole-like quasiparticles are non-trivial. We consider 
here the case that the segment / contains no disorder, 
but a sufficiently steep tunnel barrier which makes no 
difference in the tunnelling probability of electrons and 
holes. In this case we can omit the energy dependence 
in the matrices ri.2 and ti^2- In what follows we take 
advantage of the polar decomposition 





«3 




f Ul 






r2j 









u\ 
VI 



(34) 

where m/, vi are some unitary matrices, which depend 
on a particular realisation of the barrier, and F is the 
diagonal matrix of the tunnelling probabilities F^ . Time- 
reversal symmetry in the segment / is assumed. Once 
the dependence on energy in the matrices m/, vi and 
F is disregarded we obtain from Eqs. (|33| , |3^ ) the right 
refiection matrix 



TR = 



ui 
u 



COS! 



ie^^ sin 6 



uj 
u\ 



, (35a) 



'•I / \ '^"^ " ^ " / \ " , 
sin 61 = F [(l-e2*'^-*(l-F))(l-e~2*"^^(l-F))]"^ (35b) 



(1 



)[1 



■(1-r)]-^ 



(35c) 



The left matrix is given by Eq. (|5|) and describes 
the refiection from the disordered wire. Taking advantage 
of the polar decomposition we can write 



uo 
u* 



e"*' 




pi4> 



< 
ul 



(36) 



where uo is a random unitary matrix and is the di- 
agonal matrix of the eigenphases. We see that all in- 
formation contained in ui disappears statistically from 
the eigenvalues of r^rfl because the product uJmq can 
be regarded again as a random unitary matrix. Thus 
the disorder-averaged LDOS depends only on the trans- 
mission eigenvalues F^- of the tunnel barrier. Below we 
calculate the mean LDOS for a single-channel wire and 
for a multi-channel wire provided the tunnelling proba- 
bilities are the same for all channels, i.e., F,- = F. 



B. Single channel wire 

We start with the calculation of the mean LDOS for 
e < A in the case of the single- channel wire = 1. For 
e < A the phases x and 6 defined in Eq. ( psf ) are real and 
both tl and rn are unitary 2x2 matrices. We denote 
u^uq = exp(i'0), where tj) \s a, random phase distributed 
uniformly in the inter val ( 0, 27r). We insert the refiection 
matrices from Eqs. (35a, 3^) directly to Eq. (p^). The 
matrix (1 — r^rn) can be easily inverted. Taking the real 
part we notice that the zeroes of det(l — r^rii) define 
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FIG. 3: The mean LDOS (^) in a single-channel disordered 
wire at ballistic distance from an NS interface of finite trans- 
parency. The parameters Ats and F = 0.5 are chosen to fix 
the combination Ats{2 — F)/F « 1.5. The curves are close to 
each other for e ^ A where Eq. (M) is applicable. 



the exact positions of the quasiparticle bound states for 
e < A. The result reads 



n{x, s) 



x) S{cos9 costp 



X)), (37) 



where the argument of the Dirac (5-function corresponds 
to the quantisation condition for the bound states. The 
mean LDOS is given by the average over the phase (f> 
with the probability density of Eq. (]l9|), and over 

the uniformly distributed phase ip. The integration over 
^p is readily done with the result 



TT-e-x 



n{x, e) 



d(f) pe{4)) 



sin(0 + x) 



'-i^ + x) 



(38) 



In the limit F ^ 1 of the vanishing tunnel barrier one 
observes that x ~^ 7r/2— 0^ and — + 7r/2, so that the area 
of the integration in Eq. (|3^) shrinks to the small vicinity 
oi — (j)A and the function Pe{(t>) can be substituted by 
its value in this point. The integral approaches tt and we 
recover the result of Eq. for the ideally transmitting 
interface. 

In the opposite extreme of a high tunnel barrier (F — > 
0) both 6 and x go to zero, so that the integration area is 
not restricted and the value of the integral tends to unity 
because of the normalisation condition for the probability 
density Pe[(j))- In the limit e <C A we can set x = ^-nd 
reduce Eq. (Bq) to the following form 



n{x, e) 



:Re J dte 





sin^ 9 



■t{t- 2i£Ts) 



(39) 



where sinfl ~ F/(2 — F), according to Eq. (35b). From 
Eq. (^) we find that 



n(x^ e) = TTETf 



2-F 
2F ' 



ST. < 



F 



2-F 



(40) 



which coincides for F = 1 with the result of Eq. ( pOj ) for 
TV = 1 . In the dirty limit A ^ and for a high tunnel 
barrier F ^ 1 the result of Eq. ( p9| ) is applicable almost 
up to the value of e = A. It describes the formation of 
the pseudo-gap near the energy t~^F due to the normal 
reflection from the barrier. 

The exact expression (^) additionally accounts for the 
peak at e ~ A. This expression can be further gener- 
alised for energies higher than A by means of the analyt- 
ical continuation e = iuj, with the result 



n{x,e) = — Rc / da pi{(t) 



smhQ{a) 







sinli^ Q{a) + siii^ e 



1 a^juj) + 1 - F 

Q Cr) = o 1 ^ 2/ \(1 FT 

2 1 + a'^[uj)[l — I ) 



im- 

2 a 



(41a) 



(41b) 



where the function pi(cr) — (1/2) exp(— (t/2) is the prob- 
ability density (^4|) for a single-channel wire, the function 
a{uj) is deflned in Eq. (21), and the continuation to the 
real energies uj —^ —ie + 0'^ is performed. 



C. Multi-channel wire 

The disorder-averaged LDOS for ^ 1 can be found 
straightforwardly for the case of equivalent tunnelling 
probabiliti es F, = F. Then the diagonal matrices 6 and 
X in Eqs. (35b, 35c) can be regarded as scalars. It is con- 
venient to make use of the analytical continuation e — iuj 
and define 



a{Lu) 



P - 



1-f a2(,^)(l_F)' 



(42a) 
(42b) 



where R — diag {Ri, ... Rn) is the diagonal matrix of 
refiection probabilities for the disordered wire with a fic- 
titious absorption uj. In the parametrisation (22) the 
joint probability density of Rj is related to the orthogo- 
nal Laguerre ensemble (|2^). Note that the quantities p, 
Rj and a(w) take real values in the interval (0, 1) when 
UJ is real. 

The basic expression (|l|) for the mean LDOS is man- 
ifestly invariant under an arbitrary unitary rotation of 
the matrix product r^rfl. From Eqs. (p5a,pq) we obtain 



fpR 




U = UqUjuJuo, 



V cos6' [/ -isin6' 
)\ -isinO cos9 W 

Uq = diag(uo,tiS), 



, (43a) 
(43b) 



where we take advantage of the quantities defined in Eq. 
( ^ ) . The matrix uq is a random unitary matrix which is 
uniformly distributed in the unitary group (provided the 
weak disorder kpt^ 1). Hence by construction ( |43bD , U 
is th e un itary symmetric random matrix. We substitute 
Eq. (43a) into Eq. (Ilj) to express the mean LDOS as 



n{x, e) — — Re Tr 
^ ' N 



1-pR 
l+pR 



F(cos ( 



(44) 
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where 



F{z) = 



l-z(. 
pR 



Co = 



1 



(45a) 



(45b) 



\+pR 1+pR 

The average over disorder in Eq. (^J) is decoupled into 
two independent steps: the average (. . .)[/ over the group 
spanned by the unitary symmetric matrices and the aver- 
age (. . .)h over the orthogonal Laguerre ensemble of the 
reflection eigenvalues Rj. 

In the case of the finite number of channels the calcula- 
tion of average over the unitary matrices U is technically 
difficult and cannot be done analytically. However, for 
the diffusive wire, iV ^ 1, the calculation can be done 
by means of the diagrammatic technique developed in 
Ref. 

Let us briefly quote the basic substitution rules of the 
diagrammatic technique 



U*. 



(46) 



Here the matrix element Uij is represented by the black 
and white dot connected by the dashed line. The black 
dot stays for the first index i and the white dot for the 
second index j. The conjugated matrix U* is marked by 
a star. The other matrices are denoted by thick solid 
arrows. The summation over a matrix index in a dot is 
indicated by the attachment of a solid line. The aver- 
age over the unitary symmetric matrices is symbolically 
performed by pairing in all possible ways all black and 
white dots belonging to U to all black and white dots 
belonging to U* . This pairing is denoted by the thin 
solid line, which corresponds to the Kronecker symbol. 
The result of the averaging is found by inspection of the 
closed circuits in the diagram which consist of alternat- 
ing thick and thin solid lines (T-circles). Each diagram 
is weighted by a factor, which is obtained by inspection 
of the closed circuits of alternating thin solid and dashed 
lines (U-circles) . 

We expand the matrix F{z) ( p5| ) into a geometric se- 
ries and keep only the terms with equal number of U 
and W matrices. In the large-iV limit we have to take 
into account the diagrams with the largest number of 
T-circle£]. This amounts to the summation of the 'rain- 
bow' diagrams, or diffusion ladders, depicted symboli- 
cally in Fig. IJ. The corresponding Dyson equation is 

{F)u ^ i + z^iCi{F)u + zi:2C2{F)u, (47a) 

oo 

El = ^ Wnz"'' [Ti-C2{F)ur [TrCi(F)y]"^i , (47b) 

oo 

S2 = 5] W^„z"-i [TTCi{F)ur [TTC2{F)ur-' , (47c) 



where the weight factors 

W„ = jvi-^"(_i)"-i + o (iv-2n^ 

n](n — 1)! 



(48) 




= W, + z W2 



: W, 



E2 ^/C,<F>,/C, V^,<F>VC, VC,<F>V^2 V^,<F>VC, 

FIG. 4: Diagrammatic representation of the Dyson equation 
(0) for {F{z))u. 



have been found in Ref. |3^. Taking the coefficients Wn to 
the leading order in N we define the generating function 



h{s) 



E 

n=l 



(49) 



which may be used to reduce Eq. (^ ) to 

(F(z))^ = i + z2/i(z2siS2) (S1C1 + S2C2) (F(z))^, 
si,2 - TrCi,2(F(z))^. (50) 

The matrix {F(z))u has to be eliminated from the 
Dyson equations (|50|). After that it is very convenient 
to transform to the new scalar variables 



X 



S2 " 51 

N 



Y 



which obey the equations 
X + 1 1 



with 



2 

sin^ 



f{X,Y) 



-Tr- 



52 + Si 

N ' 



1 



N" l+pRfiX,Y) 

x'^cos^e = 1, 

{l-X){Y + X) 
{1 + X){Y-X)' 



(51) 

(52a) 
(52b) 

(53) 



where we hav e sub stituted z — cos 9 and the matrices Ci , 
C2 from Eq. ( |45b| ). 

In terms of the variables X and Y the mean LDOS 
is simplified to 



n{x,e) — Re X{uj) 



(54) 



where the bar stands for the average over the ensemble 
of the reflection probabilities: X = {X)i^. 

Let us first consider the case of equal reflection prob- 



abilities Rj = R. The matrix U in Eq. ( |45D commutes 
with Ci and C2 and can be diagonalised, hence the prob- 
lem becomes equivalent to that of a single channel wire. 
The solution of the self-consistent equations Eq. (^2|) is 
given by 



X = - 



sinh Q 



Vsinh^ Q + sin^ ' 



= -lnpi?, (55) 
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FIG. 5: The mean LDOS in a diffusive normal-metal wire 
in the vicinity of an NS interface of finite transparency. The 
curves are calculated from Eqs. ( ^4|j57| ) for Ats = 1 and the 
effective tunnelling probability F varying from 0.1 to 1 in steps 
of 0.1. The thick line corresponds to F = 0.5. 



which coincides with the result of the exact integration 
over U. This proves that the set of diagrams which we 



took into account in Eqs^47|) is sufficiently complete. 
The substitution of Eq. (^) in Eq. (54) and the addi- 
tional average over the reflection probability of a single 
channel wire yields the mean LDOS of Eq. (|4l|). 

In the multichannel (diffusive) limit ^ 1 the re- 
flection probabilities Rj are, in fact, not equal. More- 
over they effectively repel each other according to Eqs. 
( p^ , p3[ ). In this case Eq. ( |5^ ) can no longer be solved 
in closed form. In other words, the averages over the 
random matrix U and over the reflection eigenvalues Rj 
cannot be performed separately. 

In order to proceed one has to take advantage of the 
self-averaging property of the variables X and Y in the 
limit N ^ 1. Indeed both variables are deflned via the 
traces si.2 and can be thought as the arithmetic means 
of N fluctuating quantities. From physical point of view 
the variable X is proportional to the one-point Green 
function, therefore it is self-averaging in a diffusive metal. 

Thus we can construct the self-consistent equation for 
X by taking the average over R on both sides of Eq. 
(HI). We assume a fixed value of f{X,Y{X)) = f(X) 
on the right side, neglecting the fluctuations oiX. Taking 
advantage of the square-root approximation ( [26| ) of the 
density p{a) we obtain 



2ujTs + C 



2 "STry VC 2c.r.+ (l+p/(X))c' 



(56) 
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FIG. 6: The mean LDOS in a normal-metal wire in the 
vicinity of an NS interface of finite transparency. The dot- 
ted curve is found from Eq. ( ^l| ) for the single-channel wire 
A'^ = l. The solid curve is calculated from Eqs. (^,^^ for the 
diffusive wire A'^^1. The dashed curve represents the result 
of the Usadel equation, Eqs. ( ^5|j69| ) , calculated for the cor- 
responding value of the parameter 73 — Ars(2 — F)^/(2F)^. 
The figures show the energy dependence of the mean LDOS 
for the dirty Ats = 0.3 (top) and the clean Ats = 5 regime. 



explicitly giving rise to the equation 
{a'{uj) + l-T) {Y 



{l + a^[u){l-T)){Y{X)-X) 
2 



1 



l + X 




\+X + UJT. 



(57) 



The integral on the right-hand side can be carried out 



which is an algebraic equation for X. It can be analyti- 
cally continued to real energies oj — —ie -I- 0+ and solved 
numerically by iteration. The disorder-avera ged LDOS 
is determined, then, from Eq. (^). Equation (p^) is ob- 
tained in the quasiclassical limit of a large number of 
channels. This result does not change if we neglect that 
U is symmetric or take the unitary Laguerre ensemble in 
Eq. (p3| ) instead of the orthogonal one. 

The weak-localisation correction (which we simply de- 
fine as 1/iV correction) can, in principle, be determined 
within the present approach. It has three different 
sources. First of all an additional class of diagrams, 
namely the cooperon-like diagrams, have to be taken into 
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account in the Dyson equation (^) . Secondly the term of 
sub-leading order in the large- expansion of the weight 
factors Wn has to be included. Finally the correction of 
order 0{N~^) to the limiting form ( pq ) of the probabil- 
ity density p{a) has to be considered. The calculation of 
the weak localisation correction to the LDOS is, however, 
beyond the scope of this paper. 

In some limiting cases Eq. (^7|) allows for a transparent 
analytical solution. In the absence of the tunnel barrier 
r = 1 we obtain 



2 

X = 1- 2 " ^ (i + n- VnV^Tn 

1 -I- a-^ V 



, (58) 



which coincides upon the substitution in (54) with the 
result of Eq. (|25| ) in the large-N hmit. 

In the hmit Ar^ < T^, E^^^ (|5^) leads to the BCS result 
for the local density, Eq. (|27|). 

For small energies, w <C A, we can put a{uj) — 1 and 
obtain 



X 



I sm 



2sin^ 



The mean LDOS (M) for e < A is then given by 



n{x, e) — Re 



sm 



4 sin 



(59) 



(60) 



with sin( 



-r/(2-r). 

\-2 



This result describes the scaling 



Eg ~ T~^r^(2 — F)^"^ of the size of the pseudogap Eg 
with the transparency of the tunnel barrier F, which is 
illustrated in Fig. |^. We observe that in the limit F^ ^ 
Ats <C 1 two different types of bound states contribute to 
the LDOS at energies below A. One group of the bound 
states is responsible for the monotonous increase of the 
LDOS to its bulk value at the scale t~^F^ while another 
group gives rise to the formation of the peak near e = A. 



V. USADEL EQUATION 

The aim of this Section is to compare our results in the 
limit ^ 1 to the results of the conventional quasiclassi- 
cal theory based on the Usadel equation. It is important 
to remember that the Usadel description is justified only 
in the dirty limit Ats <C 1, while it is not restricted to 
the clean superconducting material as is the case with 
our calculation. In the quasiclassical context the super- 
conductor as well as the normal metal are characterised 
by their diffusion constants Dg, D„ and normal-state re- 
sistivities ps, Pn, which are combined into the mismatch 
parameter 



7 ^ 



Pn^n 



(61) 



where ^n,s = \/Dn,s/^ are the diffusive coherence 
lengths. Hence, the comparison has to be done in the 
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FIG. 7: The mean LDOS in vicinity of an NS interface 
of finite transparency calculated from Eqs. ( ^^Js^ ). The pa- 
rameters Ats and r = 0.1 are chosen to fix the combination 
7I = ATsi2 - rf/(2r? ^ 1.06. The curves coincide for 
£ ^ A, where Eqs. ( pOUtiI) are applicable. The result of the 
Usadel equation, Eqs. ( ]6q| , p9| ), is indistinguishable from the 
dashed line. 



limit 7^1, where the "rigid" boundary condition is 
valid. 

In the case of the perfectly transparent NS bound- 
ary and vanishing mismatch parameter the_LDOS at the 
interface found from the Usadel equatior£3li3 is simply 
given by the standard BCS formula and, therefore, coin- 
cides with our expression ( p7| ) in the dirty limit At^ <C 1. 
Thus, there is not too much to compare for the case of 
transparent boundary. However, if the NS interface is 
not perfectly transparent (F < 1), even the limit of small 
mismatch parameter 7 is not completely trivial. Let us 
now discuss the Usadel equation for this case in some- 
what more detail following the calculation of Ref. [ij. 

The transparency of the interface enters the theory 
through the parameter 



7s 



Rb 

Pn^n 



(62) 



where Rb is the product of the barrier resistance and its 
area. The Usadel equation in the normal metal {x < 0) 
takes the form 



e::(a;)-tjsine„(a;) -0, 



(63) 



where uj — —ie + 0~^ is the imaginary energy, while in the 
superconductor (x > 0) the equation reads 

^e';(a;) -wsine,,(a;) + A(a;)cose,(a;) = 0, (64) 

where A(a;) is the gap function. (For a sake of simplicity 
we restrict ourselves to zero temperature.) The func- 
tions G(x,x) = cosO„_s(a;) and F{x,x) = sin8„^s(x) 
parametrise normal and anomalous quasiclassical Green 
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FIG. 8: The mean LDOS from the random matrix theory 
(solid lines), Eqs. ( p^j57| , is compared to that from the Us- 
adel theory (dashed lines), Eqs. (^, ^), for the corresponding 
value of the parameter 7I = Ar42 - V'^)l{2Vf. The curves 
always coincide for small energies e <C A. The perfect agree- 
ment in the entire energy range is found in the dirty limit 
Ats 0, where the Usadel equation is justified. 



functions in energy representation, averaged over angle 
and disorder. The LDOS near the interface is given by 



n(0,e) = Recose„(0). 



(65) 



Far away from the NS interface the Green functions 
aquire their bulk values 



cos0„(— 00) = 1, 



cos 0s (00) 



VA2 



(66) 



The finite transparency of the interface comes ipto play 
in the appropriate matching conditions at x = 0ll3 

TsCneUO) - sin [9^(0) -e„(0)], (67a) 
liu^'M - e.©;(0). (67b) 

Once the superconductor is sufficiently clean the first 
term in Eq. ( p^ ) can be disregarded, hence Os(x) 
Os(oo) and A(a;) = A for a; > 0. This justifies the "rigid" 
boundary conditions, which are used throughout the ar- 
ticle. 

The first integral of Eq. (^3|) is readily found 

Dn. 



^(e;(x))2+c^cose„(a;) 



const, (68) 
where the constant is determined from the condition at 



X = — cjCL.and equals to. With the help of Eq. (67a) one 
obtainsEl 



sin2(e„(0)-es(0)) , uj 



47I 



+ -(cose„(0)-l)-0, (69) 



where 8s (0) is substituted by Gs(oo) due to the "rigid" 
boundary condition. In the limit w <C A the equation is 
simplified to 

f cos e„(0) -^Y = ^ (1 - cos e„(0)) . (70) 



n(0, e) = Re 



1 



(71) 



which is manifestly equivalent to Eq. ( |60| ) and establishes 
the following relation between the parameters 



7I = Ar, 



(2-r)^ 
4r2 



(72) 



This relation also follows directly from the definition of 
7b, up to a numerical factor, since one can effectively 
substitute Bb = (^Ve)(2 - r)/2r, p„ {h^/e)£-'^, and 

We conclude that the LDOS obtained from the Usadel 
equation always coincides with that found from Eq. ( [S^ ) 
for small energies £ ^ A. We also demonstrate numer- 
ically in Figs. ^, 0, and || that our result for ^ 1 is 
perfectly consistent with the Usadel theory in the dirty 
limit Ats <ti 1, where the latter is justified. 

One should note, however, that the agreement with the 
quasiclassical theory becomes better with the increasing 
barrier height. Indeed, in the perfectly transparent in- 
terface r = 1, the agreement is reached only in the ex- 
tremely dirty limit Ats 0, while for smaller values of 
r the dirty-limit condition is less restrictive (see Fig. |^). 



VI. CONCLUSION 

In conclusion, we compute the mean LDOS in a 
normal-metal disordered wire in the immediate vicinity 
of an NS interface at zero temperature and zero magnetic 
field. Our calculation is based on the scattering approach 
and takes into account the spatial phase coherence in the 
normal-metal. 

We derive the general formula (|l^), which expresses 
the one-point Green function in terms of the reflection 
matrices. The formula can be applied in order to calcu- 
late the LDOS (and its distribution) in the wire at arbi- 
trary distance to the NS interface. In this paper we only 
consider the mean LDOS at the ballistic distance to the 
interface so that it does not acquire a spatial dependence. 

We obtain the relation (|l|) between the disorder- 
averaged LDOS near the ideal NS interface and the prob- 
ability density of the eigenphases of the matrix correlator 
'''o{£)fo{~sy , where ro{e) is the reflection matrix for the 
semi-infinite normal-metal wire. 

We also study in detail the case of the normal- 
superconductor tunnel junction and derive the self- 
consistent equation ( ^7|) that determines the LDOS in 
the diffusive normal metal. In the dirty limit our ex- 
pression coincides with the LDOS found by Golubov and 
Kupriyanovlla from the Usadel equation. 

The quasiclassical analysis of the Green function at the 
NS interface o f finite , teassparency has been performed by 
many author alSMH'LBcil in connection with the bound- 
ary conditions of semiclassical superconductivity. How- 
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ever, to our best knowledge there exists no counterpart 
to the equation ( p7|) in the hterature. 

In the case of an ideal NS interface the LDOS is found 
to be almost independent of the number of channels in 
a wire, except for very small energies, hence its insen- 
sitivity to phase-coherence effects. This persists to the 
case of finite transparency provided the clean limit con- 
dition Ats 3> 1. In the dirty limit Ats <C 1 and small 
transparency F ^ 1/A^ the situation is different and the 
phase-coherent effects play a role. 

The effect of Anderson localisation is seen in the linear 
increase of the LDOS, n = (7r/4)(iV 1)£Ts(2 - r)/r, for 
energies lower than Ec = l/N'^Tg. In the diffusive metal, 
oo, the LDOS increases as the square root of energy 



n = Re y/—ieTs{2 — r)/r. The form of the crossover in 
energy dependence of the LDOS from linear to square- 
root behavior is given by Eq. (^5|) for weak disorder and 
perfect NS interface. 
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